No naked singularities in homogeneous, spherically symmetric bubble spacetimes? 
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We study the evolution of bubble spacetimes in vacuum and electrovac scenarios by numerical 
means. We find strong evidence against the formation of naked singularities in (i) scenarios with 
negative masses displaying initially collapsing conditions and (ii) scenarios with negative masses 
displaying initially expanding conditions, previously reported to give rise to such singularities. Ad- 
ditionally, we show that the presence of strong gauge fields implies that an initially collapsing bubble 
bounces back and expands. By fine-tuning the strength of the gauge field we find that the solution 
approaches a static bubble solution. 
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I. INTRODUCTION 

Within five-dimensional Kaluza-Klein theory, negative 
energy initial data configuration at a moment of time 
symmetry have been presented in the past P, 0] ■ Among 
the reasons for considering negative energy solutions is 
that naked singularities are associated to them. There- 
fore studying such solutions is an attractive problem to 
further test the cosmic censorship conjecture jsj. Addi- 
tionally, as these bubbles can be obtained via semiclassi- 
cal tunneling (following Witten's construction it 
is important to understand their dynamical behavior. To 
obtain a better feeling for the dynamics of the solutions 
presented in which describe a bubble spacetime, Cor- 
ley and Jacobson Q analyzed the initial acceleration of 
the bubble (which can be calculated analytically with the 
information contained on a single hypersurface.) They 
found that negative mass bubbles start out expanding 
and consequently argued a naked singularity would be 
unlikely. This conjecture is backed by the calculations at 
the initial hypersurface together with the intuition gained 
from a time symmetric scenario. Namely, that if the bub- 
ble were to reverse its behavior it will go through another 
time-symmetric phase which would also suggest the 
bubble should expand (assuming the solution at this time 
of symmetry still belongs to the family presented in 0). 

Recently however, a numerical study was presented 
which computed explicitly the solution to the future of 
the initial, time-symmetric, initial data. In the nega- 
tive mass case, the solution found appeared to indicate a 
very puzzling development. Namely, that as the bubble 
expanded, it would encounter a naked singularity on its 
way. 

In the present work we reexamine the dynamical be- 
havior of the data presented in 0] via numerical simula- 
tions, and further study a more general case where the 
initial data describes an initially collapsing bubble with 
negative mass. We restrict our simulations to the con- 
text of homogeneous spherically symmetric spacetimes. 
To summarize our findings, we see that negative mass 
bubbles will not lead to a naked singularity. By choosing 
different members of the initial data family, these could 
start out either expanding or collapsing. In the former 



case, the bubble continues on expanding while in the lat- 
ter case, the bubble bounces back without ever collaps- 
ing. We additionally study the case of a positive mass 
bubble which starts out collapsing. Here, depending on 
the strength of the gauge field at the initial hypersurface, 
the bubble collapses to a black string or bounces back to 
expand "forever" (i.e. for as long as we let the simula- 
tion run). Furthermore, these two distinct possibilities 
give rise to critical behavior on the threshold of collapse 
or expansion. 



II. SET UP: BUBBLE DATA 

We consider a generalization of the time symmetric 
family of initial data presented in We start with a 
spacetime endowed with the metric 



ds^ = -dt^ + U{r)dz^ 



dr" 

W) 



(2.1) 



where dfl^ = dd^ + sin^ ddip"^ is the standard metric on 
the unit two-sphere and U{r) is a smooth function 
that has a regular root at some r = > 0, is everywhere 
positive for r > r_|_ and converges to one as r — > oo. The 
coordinate z parametrizes the extra dimension which 
has the period P ~ 47r/C/'(r+). The resulting spacetime 
{t,z,r > r_|_,z?, (^} constitutes a regular manifold with 



the topology 



X S . The bubble is located where 



the circumference of the extra dimension shrinks to zero, 
that is, at r = r4.. 

Additionally, we consider the presence of a U(l) gauge 
field in the following form. 



Aadx" = fc(r+" - r-")dz. 



(2.2) 



with k an arbitrary constant and n an integer greater 
than one. It is important to point out two things. First, 
the gauge field here considered has a more generic depen- 
dence than those previously considered in ^ , where only 
the case n — 2 was discussed. Second, the symmetries of 
the problem would a priori allow for the gauge field to 
have non-trivial t and r components. However, it follows 
from Maxwell's equations that those contributions would 
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lead to a Coulomb-like electric field which is divergent at 
the location of the bubble. For this reason we set those 
components of Aa to zero. 

In the time-symmetric case, initial data satisfying the 
Hamiltonian constraint obeys 



U{r) 



1 



(2.3) 



requires 



(2.4) 



where k = knj ^ {n — l)(2n — 1). The parameter m is 
related to the ADM mass via Madm = m/A. 

The fact that the bubble be located at r = 
that = U{r^) = 1 — rfi + b ~ k^, where m 
6 = k = k/r^ . We also require 

< r+U'{r+) = 2-m + 2{n- 1)P 

and avoid the conical singularity at r = r+ by fixing the 
period of z to P = 47r/J7'(r+). Repeating the analysis 
in one finds that the initial acceleration of the bubble 
area is given by 

" ~3" 



1 



-(n-l)(n-2) 



(2.5) 



For n = 2 this agrees with the formula found by Corley 
and Jacobson, and as they discussed, for negative mass 
bubbles, the initial acceleration is positive, and only if 
chosen appropriately can a positive mass bubble have 
negative initial acceleration. As we will see later, in the 
vacuum case, this acceleration remains negative leading 
to a collapse of the bubble and the formation of a black 
string. In the non-vacuum case however, the strength of 
the gauge field can modify this behavior completely. For 
weak enough k the bubble continues to collapse whereas 
when k is large the bubble area bounces back and ex- 
pands. Furthermore, by fine-tuning the parameter k we 
will show that the solution approaches a static bubble 
solution to the five-dimensional Einstein-Maxwell equa- 
tions. 

We also consider non-vacuum initial data with n > 2 
which allow for richer scenarios. A particularly interest- 
ing case is that describing an initially negative accelera- 
tion with negative mass. If this were possible, and the 
bubble continued to collapse, this would likely give rise 
to a naked singularity. Thus, we are interested in data 
for which A < 0. In view of the inequality l|2.4|l and Eq. 
(|2.5|) this means that 

2 + 2(n - l)fc2 > m > 1 (71 - l)(7i - 2). (2.6) 

3 

From this we see that for 71 = 2 the mass can only be 
positive. In contrast to this, for 71 > 2 and k large enough 
such that 4fc^(n — l)(n — 2) > 3, negative mass bubbles 
can be considered. 



III. NUMERICAL IMPLEMENTATION 

A detailed description of the numerical techniques used 
in our code will be presented elsewhere 8] . Here we point 



out its salient features. 

1. The equations are written in first oder symmetric 
hyperbolic form, with a gauge choice that is such 
that the characteristic speeds are given by 0, 1 or 
— 1. These properties facilitate achieving a stable 
implementation of the equations. 

2. Regularity conditions are enforced at the bubble 
and employed to design a numerical scheme that 
satisfies a discrete energy estimate BT^. This 
discrete energy estimate guarantees numerical sta- 
bility for a related toy model problem. 

3. Constraint preserving boundary conditions are en- 
forced at the outer boundary to ensure no con- 
straint violating modes are fed through it 

4. The equations are implemented using second or- 
der accurate finite difference techniques. A non- 
uniform radial coordinate is employed to improve 
the resolution near the bubble. For all cases pre- 
sented here, the radial grid covered r e [1, 20] with 
a total of 8000 points. The outer boundary's lo- 
cation is chosen so as to have always less than a 
crossing time in our evolutions, guaranteeing outer 
boundary conditions (though consistent with the 
constraints) will not affect quantities at the bub- 
ble. Second order convergence is checked in all 
runs. As an example, the convergence of some of 
the simulations representing expanding bubbles is 
illustrated in Fig. |21 Similar behavior is observed 
in all caseS(19j. 



IV. RESULTS 



A. Brill-Horowitz initially expanding case: 



2) 



Here we evolve the Brill-Horowitz initial data in the 
case of vanishing gauge field. The bubble area A as a 
function of the proper time r at the bubble is shown in 
Fig. ^for different values of the mass parameter m. (Fig- 
ure |21 illustrates the convergence of the solution through- 
out the evolution). As expected, the lower the mass of 
the initial configuration, the faster the expansion. Em- 
pirically, and for the parameter ranges used in our runs, 
we found that at late times the expansion rate is given 
by 



A 
A 



r+(r = 0) 



(4.1) 



where a dot denotes the derivative with respect to proper 
time T. In particular this approximation is valid for the 
bubble solution exhibited by Witten P| which describes 
the time evolution in the case tti = 0. We monitored 
several curvature invariant quantities (the Kretschmann 
invariant in particular) for our numerical spacetimes and 
found no sign of divergence. 
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FIG. 1: Bubble area vs. proper time at the bubble. In 
this and the following plots, we set r+ = 1 (by choosing b 
appropriately). The figure shows four illustrative examples of 
bubbles whose initial acceleration is positive. As it is evident, 
the expansion of the bubble continues and the difference is 
the rate of the exponential expansion. The relative error in 
these curves is estimated to be well below 0.001%. 



mass is positive, one expects this singularity to be 
hidden behind an event horizon, and one should obtain 
a black string. In fact, for the solutions which are 
initially collapsing and which have vanishing gauge 
field, we observe the formation of an apparent horizon. 
Furthermore, we compute the curvature invariant 
quantity IR\jf at the apparent horizon (as discussed 
in 17]), where / = RabcdR"''"''^ is the Kretschmann 
invariant and Rah the areal radius of the horizon. For a 
neutral black string, this invariant is 12. Figure |31 shows 
how this value is attained after the apparent horizon 
forms for representative vacuum cases (with m — 1.1 
and m — 1.99), thus providing strong evidence for the 
formation of a black string. (We also find apparent 
horizons forming, or already present in the initial data, 
when a nontrivial gauge field is considered and the 
parameters are appropriately chosen. The corresponding 
bubbles will likely collapse to charged black strings and 
will be analyzed elsewhere 0-) 



m = -0.3 
111 = -0.6 
m = -0.8 
m = 0.3 
m = 0.8 
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FIG. 2: Convergence factor (CF) vs. proper time for 
the cases illustrated in Fig. This factor, defined as 
CF = log2{\\A{4:A) - yl(2A)||/||A(2A) - A{A)\\) should 
give the value of 2 for a second order accurate implemen- 
tation. Through these runs the grid spacing is defined as 
A = 3.75 X 10~^ Clearly, the plot shows that second or- 
der convergence has been obtained throughout the runs (the 
peaks are instances with 0/0 where the expression for CF is 
ill-defined) . 



B. Brill- Horowitz initially collapsing case: (n — 2) 

Next, we analyze the Brill-Horowitz initial data for 
the case in which the bubble is initially collapsing (no- 
tice that for n = 2 this implies that the ADM mass is 
positive) . While our numerical simulations reveal that in 
the absence of the gauge field such a bubble continues to 
collapse, we also show that when the gauge field is strong 
enough, the bubble shrinks at a rate which decreases with 
time and then bounces back. 

Obviously, if the collapse trend were not halted, a 
singularity should form at the origin. Since the ADM 




FIG. 3: Rescaled Kretschmann invariant IR\jj/l2 vs. 
asymptotic time for m — 1.1 (solid line) and 1.99 (dashed 
line). The first non-zero values of the lines mark the for- 
mation of the apparent horizon. At late times, both lines 
approach the value of 1 suggesting a black string has formed. 



As mentioned, for strong enough gauge fields, the pre- 
viously described dynamics is severely affected. Figure 
21 shows the bubble area vs. proper time for different 
values of k. For large values of k the bubble "bounces" 
back and expands while for small ones the bubble col- 
lapses. There is a natural transition region separating 
these two possibilities. Tuning the value of k one can 
reveal an associated critical phenomena, the 'critical so- 
lution' being a member of the family of static solutions 
given by 




r^V{r)dn^, 



where V{r) = 1 — r_/r and U{r) = 1 — r+/r. These 
solutions can be obtained by a double analytic ex- 



pansion of the black hole solutions described by Eqs. 
(2.13) and (2.15) of Ref. ^ followed by the trans- 
formations ^ — r <—>■ r — r^, r+ r_|- — r_. 
The parameters r_ and r+ (> r_) are related to the 
period of the z coordinate and to the ADM mass via 
P = 47rr+(l - r_/r+)^/2 and Madm = r+/A. Since 
these parameters are conserved throughout evolution, 
the initial data fixes the member of the static family the 
dynamical solution approaches to. Details of the critical 
behavior exhibited by this system will be presented 
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FIG. 4: Area values vs. proper time at the bubble for different 
values of k and m = 1.1. By tuning the value of k appropri- 
ately, the amount of time that the area remains fairly constant 
can be extended for as long as desired. 



C. Collapsing negative mass case 

We here restrict to cases with negative masses that 
start out collapsing. Interestingly enough we find that 
even when starting with large initial negative accelera- 
tions, which in turn make the bubble shrink in size to 
very tiny values, it bounces back without ever collapsing 
into a naked singularity. As an example, Fig. |31 shows 
the bubble's area versus time for different values of n and 
k. The initially collapsing bubbles decrease in size in a 
noticeable way but this trend is halted and the bubbles 
bounce back and expand. Although we have not found 
a simple law as that in Eq. (|4.1|) . clearly the bubbles 
expand exponentially fast. Therefore, it seems not to 
be possible to "destroy" the bubble and create a naked 
singularity. This situation is somewhat similar to the sce- 
narios where one tries to "destroy" an extremal Reissner- 
Nordstrom black hole by attempting to drop into it a test 
particle with high charge to mass ratio. There, the elec- 
trostatic repulsion prevents the particle from entering the 
hole JJ]. 



FIG. 5: Bubble area vs. proper time at the bubble. The 
figure shows three illustrative examples of bubble with 
negative mass (m = —0.1 each) whose initial acceleration is 
negative. As it is evident, the collapse of the bubble is halted 
and the trend is completely reversed. The error in these 
curves is estimated to be well below 0.001%. 



V. CONCLUSIONS 

We have performed numerical simulations to inves- 
tigate the dynamics of spherically symmetric, homoge- 
neous bubble spacetimes. These have revealed several 
interesting features. 

In the vacuum case, negative mass bubbles do not form 
singularities, rather they continue on expanding and the 
asymptotic rate of expansion obeys the simple relation 
(|4.1|l . Positive mass cases, as pointed out previously can 
start out either expanding or contracting. In the former 
case, they continue on expanding and its asymptotic rate 
obeys the same simple relation. For the collapsing case 
we have evidences that the bubble will collapse to a black 
string. It is interesting to note here that the ratio of a 
circle at infinity over the ADM Mass is given by 



P 



167r 



Madm m{2 - m) 



(5.1) 



Since < m < 2 this value is always larger than 
that needed for the Gregory-Lafiamme "instability" jl^ ■ 
Thus, a slight ^-dependent perturbation of the scenario 
considered here could indeed give rise to a naked singular- 
ity. This would require infinite time to appear as shown 
in 16], but certainly would give rise to rich dynamics as 
illustrated in 

In the non-vacuum case, the presence of a gauge field 
gives rise to additional features. First, for the positive 
mass case, it can "halt" the collapse and subsequently 
induce an expansion. Furthermore, by fine-tuning the 
strength of the field, the solution approaches a static 
bubble. This critical solution separates the collapsing 
bubbles (which, likely, collapse to a charged black string) 
from the one that expand at late times. In a subsequent 
paper Q we will show that the dynamical behavior near 
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the static bubble exhibits critical phenomena. Note that 
clearly the solution does encounter another moment of 
time symmetry at the critical value. However, for field 
strengths larger than the critical one, this is not the case. 
The bubble bounces but the solution does not go through 
a time-symmetric phase. 

Last, the gauge field can also make negative mass bub- 
bles start out collapsing if its radial dependence is steeper 
than that originally considered in jj. Intuitively this 
happens due to the contribution of the energy in the 
gauge field overcoming the repulsive gravitational effect 
of the negative mass bubble. However, as part of the field 
is radiated away, the collapse is halted and a subsequent 
expansion results. By choosing both the strength of the 
field or its radial dependence appropriately (with n > 2), 
we could make the bubble start out collapsing at a very 
rapid rate. In these cases, the bubble shrinks to very 
small sizes but eventually bounces back without forming 
a naked singularity. 
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